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GRAVITY WAVES AND INSTABILITIES IN THE LOWER AND 
MIDDLE ATMOSPHERE 

Jurgen Klostermeyer 

Max-Planck-Institut fur Aeronomie, 3411 Katlenburg-Lindau, FRG 


1. INTRODUCTION 

The atmosphere moves ceaselessly on scales ranging from the dimension of the 
earth down to the mean free path of individual air molecules. An understanding 
of atmospheric motions requires the study of very specific problems that have to 
be idealized to focus attention on the basic dynamical processes, and thus are not 
faithful in detail. MST radar observations comprise large-scale processes like lon- 
gitudinally averaged wind fields, planetary and tidal waves, and synoptic weather 
disturbances, mesoscale processes like narrow jet streams, frontal zones and atmo- 
spheric gravity waves, and small-scale processes like Kelvin-Helmholtz instability and 
turbulence preferably at scales equal to half the radar wavelength. All these processes 
interact nonlinearly, e.g. winds can enhance gravity waves and gravity waves can ac- 
celerate winds (wave-mean flow interaction) or a gravity wave at a given wavelength 
and period can enhance gravity waves at other wavelengths and periods (wave-wave 
interaction). Idealizing a problem then means that some linear terms in the hy- 
drodynamic equations are neglected to suppress unwanted types of motion and that 
nonlinear terms are neglected or simplified to study particular types of interaction 
with relatively simple mathematical tools. 

It is clearly impossible to discuss the whole variety of atmospheric motions in this 
lecture, therefore we will confine ourselves to some basic aspects of mesoscale and 
small-scale gravity waves and instability mechanisms. These processes can be (and 
partly have been) investigated by single ST or MST radars rather than expensive 
networks of radars and other techniques that would be necessary to study details of 
large- scale processes (the physics of planetary waves, for example, depends strongly 
on the zonal flow averaged around a latitude circle and thus cannot be investigated 
by measuring only the time- averaged flow over a single radar station). 

Internal gravity waves with wavelengths between ten and less than one kilome- 
ter and periods between several hours and several minutes appear to play a central 
role in atmospheric wavenumber and frequency spectra (VANZANDT, 1982). There- 
fore we discuss the propagation of gravity waves in simplified atmospheric models 
in Section 2. Their interaction with the wind, their mutual interaction and instabil- 
ity mechanisms based upon these processes will be treated in Sections 3 and 4. All 
sections stress the theoretical aspects but are completed by MST radar observations 
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showing the relevant hydrodynamic processes. Such a short treatment necessarily 
must be incomplete and sketchy but hopefully will stimulate further studies of dy- 
namic processes by means of MST radars and other methods. 

2. ATMOSPHERIC GRAVITY WAVES 


To study atmospheric gravity waves it is convenient to consider an unbounded 
non-rotating model atmosphere without molecular viscosity and thermal conductiv- 
ity. Thereby we eliminate unwanted types of dynamic processes like gyroscopic or 
inertial waves, viscosity waves and heat conduction waves resulting, respectively, from 
the Coriolis force and the dissipative terms in the equations of momentum and energy 
which play no role at the wavelengths and frequencies considered below. For fur- 
ther details of these wavetypes see VOLLAND (1969) and LEBLOND and MYSAK 
(1978). The remaining waves that can propagate through the model atmosphere 
fall into two categories: High-frequency acoustic waves due to the compressibility of 
the air and low-frequency internal gravity waves due to the gravitational force. The 
periods and wavelengths of acoustic waves are smaller than about 1 s and 300 m, 
respectively, whereas the corresponding values for internal gravity waves are larger 
than about 300 s and 300 m. The dispersion curves of both wave types thus are 
widely separated in the frequency- wavenumber plane and interactions can in general 
be neglected. We can therefore eliminate acoustic waves by assuming that the model 
atmosphere is incompressible so that the sound speed is infinite and the density 
cannot vary along the path of an air parcel: 



(i) 


where d/dt = d/dt + u- V and p and u denote density and fluid velocity. It should be 
mentioned, however, that this approximation also filters out long-period Lamb waves 
which propagate horizontally at the speed of sound and sometimes play a dominant 
role in the atmospheric response to tropospheric excitations (KLOSTERMEYER, 
1977; LIU et ah, 1982). 


Together with (1), the continuity equation and the Euler equation yield a closed 
set of equations for describing the dynamics of the model atmosphere: 


V • u = 0 

(2) 

du _ 

(3) 

P ~dt + V P “ /’g = 0 


where p and g denote pressure and gravitational acceleration. 

Often it is useful to compare the actual atmospheric state with a reference state 
defined by hydrostatic equilibrium 
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p = Po(*), P = Po(*), u = 0 ( 4 ) 

satisfying (1) and (2) while (3) reduces to 

Po* = ~Po9- ( 5 ) 

Here and in the following, the subscripts x, z and t denote partial derivatives with 
respect to horizontal and vertical coordinates and time, respectively. Equation (5) 
yields p 0 (z) for any given density distribution po{z). An important quantity of the 
equilibrium state is the Vaisala-Brunt frequency N given by 

N 2 = -g^. ( 6 ) 

Po 

It will be shown below that N is the angular frequency of an air parcel when it is ver- 
tically displaced from its equilibrium position to a slightly higher or lower level. The 
period 2n/N is about 5 to 10 min in the earth’s atmosphere. N 2 < 0 (convection) 
means that the air parcel would continue to rise or fall rather than oscillate around 
its equilibrium position. Such unstable situation in general cannot persist because 
this so-called static instability is eliminated by strong vertical mixing as fast as it 
forms. 

Small departures from the basic state (4) can be described by the perturbation 
density pi, pressure p\ and velocity Ui defined by 

p — Po + pu P — Po + Pi? u = Ui. (7) 

Substituting into (l)-(3) and neglecting terms containing productes of perturbation 
quantities yields 


Pit + tflPOz = 0 

(8) 

« 

c 

1! 

o 

0) 


u if + — Vp! - —g = 0 (10) 

Po Po 

where w-\ is the vertical component of Uj. Equation (10) indicates that a reduced 
gravitational acceleration (or buoyancy) gpi/po and a modified pressure p\ act upon 
a displaced parcel. In a stably stratified atmosphere, buoyancy pulls back any parcel 
to its equilibrium position giving rise to internal wave oscillations. Linearization is 
primarily an approximation dictated by our inability to treat nonlinear problems with 
adequate mathematical tools. It is valid for wave motions of infinitesimal amplitude 
implying that the fluid velocity must be much smaller than the phase velocity, and 
clearly filters out all nonlinear interactions between waves of different wavelengths 
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and periods as well as self-interaction. 

It is convenient to eliminate the variables p\,p\ and u x (the horizontal component 
of Uj) yielding an equation for uq alone: 

N 2 

+ JV^Vjuq = 0 (11) 

9 

where V h = (d/dx, d/dy, 0). For constant Vaisala-Brunt frequency N, (11) has 
plane wave solutions of the form 

u>, = Wie i{lx+M *- w (12) 

Then the dispersion relation 

+ (U) 

must hold for a nontrivial solution. 

The imaginary term in (13) is a consequence of the fact that for constant N, the 
unperturbed density p 0 decreases exponentially with height while the wave energy 
density should remain uniform in the absence of energy sources. Equation (6) in fact 
yields 

Po{z) = Po(0)e~% (14) 

where H = g/N 2 is the density scale height. To derive the perturbation energy 
density we note that the vertical displacement ( of a fluid parcel from its equilibrium 
level is given by 

Ci = (15) 

so that integrating (8) with respect to time yields 

Pi = (16) 

From the scalar product of Uj with (10), and from (15) and (16) we then obtain 

’?(U? + N\ 2 )} +V-(p 1 u 1 ) = 0 (17) 

l Z j ( 

where po\x\/2 and p 0 N 2 ( 2 / 2 are the kinetic and potential wave energy densities and 
puj is the energy flux up to an arbitrary nondivergent contribution. The average 
kinetic and potential energy densities are equal, and the average wave energy density 
E varies as (ECKART, 1960) 

E oc p 0 (z)Ni| 2 - (18) 

Constant E thus requires 
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wi oc e 3 ^. (19) 

Indeed, with M = ReM -f ilmM we obtain from (13) 

Im M = (20) 

and with m = ReM 


u > 2 = 


/ 2 + m 2 + 


JV 2 


( 21 ) 


The frequency is zero for vertical phase propagation (/ = 0) and equal to N for 
horizontal phase propagation (m = 0) if l 2 » 1/4 H 2 . For cj > N y no waves are 
possible. The hyperbolae (dashed curves) in Figure la represent contours of constant 
u; in the wavenumber domain. The group velocity 


duj dus 

Cj = ( aT’ 


( 22 ) 


is normal to the contours always pointing away from the ordinate. It is typical of 
internal gravity waves that the vertical components of the phase and group velocities 
have different signs. 



Fig. 1. (a) Contours of constant frequency in the wavenum- 
ber domain for internal gravity waves in an incompressible 
fluid (dashed lines) and in a Boussinesq fluid (continuous 
lines). The group velocity c g is normal to the contours, (b) 
Internal gravity waves are transverse with the fluid velocity 
parallel to the wave fronts. 

For short wavelengths with l 2 + m 2 >> 1/4JT 2 , (21) yields 


J 1 - N 2 cos 2 0 


(23) 
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where without loss of generality, 6 is the angle between the horizontal and the 
wavenumber vector k = (/,m) with }9\ < 7r/2. The short wave dispersion equation 
(23) could have been obtained directly from the equations of motion by applying the 
Boussinesq approximation which consists of neglecting density variations in the in- 
ertia terms but retaining them in the buoyancy term. Then the last term on the left 
of (11) vanishes and (23) follows immediately. With the Boussinesq approximation, 
the hyperbolic contours u) — const in Figure la degenerate to linear asymptots of the 
hyperbolae including the angle 0 with the l axis indicating that this approximation 
is the better the larger k. 

From (9) and (10) we find that plane internal gravity waves are transverse with 
Ui in the x — z plane parallel to the wave fronts (Figure lb). Note that at 6 = 0, 
the phase propagation is horizontal and U! vertical, showing that for small vertical 
displacements and N 2 > 0, an air parcel oscillates at the Vaisala-Brunt frequency. 

A concise description of further details of the linearized theory of wave propa- 
gation such as wave energy, wave action, slowly varying wave trains in nonuniform 
media etc. is given by BRETHERTON (1971). 

We will complete this chapter by showing MST radar observations that indicate the 
relevance and applicability of the foregoing to dynamic processes in the troposphere, 
stratosphere and mesosphere. Figure 2 contains contour lines of the power density 
of vertical velocities observed after the passage of a severe thunderstorm. Convective 
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Fig. 2. Time-height contours of power spectra of vertical ve- 
locity measured with the SOUSY VHP Radar on 2 June 1978 
after the passage of a thunderstorm. The peaks of the spec- 
trogram correspond to a power of 1.1 x 10 _5 m 2 s~ 2 . The dot- 
ted curve represents the height profile of the Vaisala-Brunt 
period obtained from radiosonde data (from ROTTGER, 
1980). 
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activity in the troposphere at times of thunderstorms should in particular excite inter- 
nal wave motions in the stratosphere because rising columns of tropospheric air may 
penetrate a short distance through the tropopause and transfer kinetic energy to the 
stably stratified stratosphere (STULL, 1976). The spectrogram in Figure 2 demon- 
strates that oscillations with periods larger than 2 tt/N are observed at all heights 
with the strongest oscillations occurring above the tropopause at a height near 10 
km. No oscillations are found at periods less than 2 n/N so that it appears possible 
to derive mean potential temperature profiles from the cut-off period of gravity wave 
spectra (ROTTGER, 1980). 



Fig. 3. Time variation of vertical velocity at four heights 
measured with the mobile SOUSY VHF Radar at Andenes 
(Norway) on 21 January 1984. Spectral analysis reveals 
significant peaks at periods near 5, 16 and 30 min (from 
RUSTER, 1984). 

Figure 3 shows the time dependence of the vertical velocity component of gravity 
waves measured in the arctic mesosphere. The spectra of the time series have three 
significant peaks at periods near 5, 16 and 30 min. Spectra with multiple peaks 
appear to be the rule rather than the exception at mesospheric heights. Almost all 
radar records show the presence of horizontally propagating waves with periods near 
2n/N that must be generated by local sources rather than sources at lower atmo- 
spheric regions. Two possible source mechanism, Kelvin- Helmholtz and parametric 
instabilities, will be addressed in the next two chapters. There seem to be no detailed 
case studies of the horizontal and vertical phase propagation of short-period gravity 
waves in the literature. 

3. WAVE-MEAN FLOW INTERACTION 

In the previous chapter we discussed linear gravity waves that do not interact. 
If interaction among waves or between waves and mean flow takes place, we can 
distinguish between weak and strong interactions. In weak interactions, the space 
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and time scales of energy exchange processes greatly exceed the typical wavelengths 
and periods so that the solutions of the equations of motion can be expressed as 
power series of small nondimensional parameters and can be found by perturbation 
methods (KEVORKIAN and COLE, 1981). An example will be given in the next 
chapter. Here we will study two types of strong wave-mean flow interaction charac- 
terized by energy exchanges occurring over scales comparable to the wavelength and 
period. The first type is the critical layer absorption of internal gravity waves loosing 
their momentum to a vertically sheared mean flow, the second one is the stability of 
parallel flows to infinitesimal wave perturbations. 

3.1. GRAVITY WAVE CRITICAL LEVEL 

Besides being partially or totally reflected by a height dependent mean wind u 0 (z ), 
an internal gravity wave can meet a critical level z = z c at which u 0 (z) is equal to 
its horizontal phase speed. In a continuously stratified shear flow, the Richardson 
number is defined by 

N 2 

Ri = — (24) 

“o* 

representing the ratio of the energy required to interchange vertically adjacent fluid 
parcels against the gravity acceleration to the available kinetic energy (e.g. CHAN- 
DRASEKHAR, 1961). For Ri > 1/4, ray theory describing the propagation of waves 



Fig. 4. Schematic of progression of a gravity wave packet 
propagating toward a critical level z c at four different times. 
The horizontal component of k remains constant. The fluid 
velocity is parallel to the wave fronts. 
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in a slightly nonuniform medium can be used to study how an interned gravity wave 
packet approaches a critical level (BRETHERTON, 1966). The result is summarized 
in Figure 4 showing that the wavenumber increases and the perturbation velocity 
Ui becomes more and more horizontally oriented as the packet comes closer to z c „ 
Moreover, the time required for the packet to reach z c becomes infinite indicating 
that the packet would be effectively absorbed rather than reflected or transmitted. 

Since ray theory becomes invalid at the critical level, the Frobenius method was 
used to determine the behaviour of an internal gravity wave at z c (BOOKER and 
BRETHERTON, 1967). We again describe the perturbed state of the atmosphere 
by (7) but assume 


u = (uoCO 4- ui, ^i) (25) 

considering only two-dimensional motions in the x-z plane. Substituting into (l)-(3), 
neglecting products of perturbation quantities and eliminating P\,P\ and U\ yields 
under the Boussinesq approximation 

(P d 

— (u>i rl + U> 1 „) - — (uo„uHx) + N 2 w lxx = 0 (26) 

dV at 

where d/dt ~ d/dt + u 0 d/dx. Equation (26) reduces to the Boussinesq approximation 
of (11) for u Q = 0, and with N = const, u 0 = const and (12) yields the dispersion 
relation 


= N 2 cos 2 0 (27) 

where Q = u? - lu 0 is the intrinsic (or Doppler-shifted) frequency. If u 0 is height 
dependent, there is a new restoring effect due to the vertical derivative of the mean 
vorticity uo 7 . Assuming a solution of the form 

u>, = W x {z)e' l, < x - Ci) (28) 

then gives the Boussinesq form of the Taylor- Goldstein equation 

(ti 0 - c) 2 W? + [N 2 - (u 0 - c)u" - («o - c) 2 l 2 }W, = 0 (29) 

with the prime denoting differentiation with respect to z. The original form of the 
Taylor- Goldstein equation obtained from (l)-(3) without Boussinesq approximation 
is 


[po(uo-c)Wl]'-(pou' 0 W } y- 


+ Po(uo - c)t 2 

u 0 -c 


W x =0. 


( 30 ) 


Besides containing the important special cases of internal gravity wave propagation 
and static instability for u 0 = 0, (29) is the starting point for analysing the effects of 
wind shear on gravity waves. For further studies of the critical level problem we note 
that at z = z c ,u 0 - c = 0 so that (29) has a singularity across which solutions have 
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to be joined. We assume that around z = z c ,u 0 and N can be expanded in power 
series 


u 0 = c + u' 0 (z c )(z - z c ) + ■ ■ ■ 

N = N(z c ) + N'{z c )(z-z c ) + --- 

with Uq(z c ) ^ 0 and try to find a solution of the form 

Wi = a,{z - z c ) 0 ‘ +; , a 0 ± 0. 

j=o 

Substituting (31) and (32) into (29) yields 

^ ± t 1 ~ («(*.) “ j)* 

requiring Ri{z c ) > 1/4 for real fj.. Then the general solution near z = z c is 

Wi ss A(z - Zc)’*" 1 + B(z - 0 c)’ - " 1 . (34) 

The foregoing analysis gives no indication how to join the solutions across z = z c . 
A proper treatment requires some length (BOOKER and BRETHERTON, 1967) and 
is therefore beyond the scope of this lecture. One finds that for u' 0 (z c ) > 0, the A- 
wave W\a = A(z — z c ) l ^ 2+tti is associated with upward propagating energy for z < z c 
and z > z c whereas the B-wave W\q = B(z — z c ) 1/,2 " tM is associated with downward 
propagating energy. For u f 0 (z c ) < 0, the roles of both waves are reversed, and we 
obtain 


(31) 

(32) 

(33) 


|Wm(*>*c)| = \W lA (z < z c )\e-»*'*™°M 

\W lB (z>z c )\ = \W lB (z < z c )\e“^<^ [ i 

indicating that both waves are severely attenuated as they cross z c . The vertical flux 
of horizontal momentum is independent of height on each side of z c but discontin- 
uous across 2 C with the attenuation factor e~ 2?rM suggesting that there is significant 
horizontal momentum transfer to the mean flow at the critical level. 

From (34) we find that near z Cl 

|u»l| ~ \z - 2 C |^, |«l| ~ |z (36) 

indicating that W\ is well behaved whereas Uj tends to infinity at z c . The wave- 
associated shear tends to become infinite at an even higher rate: 

( 37 ) 

According to BOOKER and BRETHERTON (1967), these infinities are spurious 
in the sense that nature always generates wave packets of finite spectral width so 
that all Fourier components with infinitesimal amplitude encounter critical levels at 
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different heights and the infinities do not in fact occur in the total integrated distur- 
bance. Nevertheless, gravity waves generated by sufficiently monochromatic sources 
may produce large- amplitude oscillations near critical levels so that one of the hith- 
erto neglected mechanisms may become important: (1) Nonlinearities due to large 
perturbation amplitudes; (2) molecular dissipation due to strong vertical variations; 
(3) instability due to strong vertical shears. 

3.2 KELVIN-HELMHOLTZ INSTABILITY 

For illustrating some of the mechanisms and concepts of dynamic stability, we will 
work through a classic problem that was first addressed in 1868 by Helmholtz and in 
1871 by Kelvin, and demands only little mathematics. We consider a basic flow of an 
incompressible inviscid two- layer fluid that has unlimited extent, constant densities pi 


z 


U 


i 


Gi 


0 


X 


— QfcsBj 


Fig. 5. Kelvin-Helmholtz flow configuration. 

and pi i, and horizontal velocities Ui and Uu (Figure 5). We assume that the flow is 
disturbed by an infinitesimal perturbation displacing the interface between the two 
layers from z = 0 to 


* = G(-M)- (38) 

Since the model is discontinuous, the Taylor- Goldstein equation (30) must be solved 
separately in each layer, and both solutions must be matched across the interface 
z = Ci by two matching conditions. 

(1) The interfacial displacement has to be continuous. To first order a'i(-fO) = 
(it + U i (lo- and iui(— 0) = ( u + UliCix at the upper and lower sides of the interface, 
respectively. For a travelling wave solution of the form (28) we then get 


a d 
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Wi(+Q) = ^i(-O) ( 39 ) 

Uj — c U 7/ -c' 

(2) At z = fj, the total pressure must be continuous. Using Bernoulli’s theorem 
(e.g. LAMB, 1945), we obtain 


Pi 


(i u , - c)u»;(+o) - Yffz~ c w ^ +0 '> 


pn 


(Un - cH(-O) - 


The perturbation is further assumed to decay to zero for \z\ 
boundary conditions 


9 

Uu-c 


u;i(-0) . 


(40) 

— ► oo yielding the 


7X5j — ♦ 0 as z — ► ±oo. (41) 

For the Kelvin- Helmholtz flow configuration (Figure 5), the Taylor- Goldstein equa- 
tion (30) reduces to 


W " - l 2 W x = 0 

and must have a solution of the form 


(42) 


W\ = 


At lz at z > 0 
Be lz at z < 0 


(43) 


to satisfy the boundary conditions (41). Equations (39) and (40) then yield two 
homogeneous linear equations for A and B which are nontrivial only if 


PiiUji ± piUi ± 
Pll + pi 


Pn ~ Pi 9_ 
Pi i + Pi I 


PiiPi 

( Pit + Pi? 


(U n - Ur? 


JL 

2 


(44) 


The first term on the right of (44) may be considered as the density-weighted mean 
velocity of the two fluid layers. The wave perturbation then travels with a velocity d 
given by 


d 2 = dl- 


PnPj 

( pn + pi y 


(U n - £/,) 2 


relative to the mean velocity where 


(45) 



is the phase velocity in the absence of any basic flow. For 

(Un - Vrf > diXLfl £ (47) 

PUP! « 

c is complex so that the wave perturbation grows and remains stationary with respect 
to the mean flow. For any given Un — Uj ^ 0, the interface is unstable for sufficiently 
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small wavelengths. If pu — pu it is unstable for all wavelengths indicating the 
stabilizing influence of a density change. Finally we note that for small density 
changes, p\ l — p] « 2 pn(pn — pi)- Then we can define a Richardson number by 


Ri _ Pit ~ Pi £ 
Pi(Uu - Ujf l 

so that Ri > - (< for stability (instability). 


(48) 


The instability of the Kelvin-Helmholtz flow configuration at high wavenumbers 
is a consequence of the fact that vorticity is concentrated at a single height. A 
physically more realistic model together with the neutral stability curve obtained by 
DRAZIN (1958) is shown in Figure 6, where 

u 0 = Au tanh(^), p 0 = e~**, Ri = ( 4 $) 

The flow is stable for Ri > 1/4. We note that the expression “Kelvin-Helmholtz 
instability” is also used for shear flow instability in models with continuous density 
and wind profiles. 




Fig. 6. (a) DRAZIN’s (1958) model, (b) 

Neutral stability curve. 

We conclude the theoretical part of this chapter by summarizing some general 
results on the stability of nonrotating stratified shear flows (for details cf. LE BLOND 
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and MYSAK (1978)). 

(1) YIH’s (1955) extension of SQUIRE’s (1933) theorem: For each unstable three- 
dimensional disturbance in a stratified fluid, there is a more unstable two-dimensional 
one propagating parallel to the flow. For this reason we have considered only two- 
dimensional waves of the form (28). 

(2) SYNGE’s (1933) theorem: A necessary condition for instability is 

\{U - Crf + d&poW)' - 2 (U - c r )poN 2 = 0 (50) 

somewhere in the field of flow. Here, c r - Rec and c, = Imc. For p 0 = const, (50) 
yields U n = 0 which is Rayleigh’s inflection point theorem. 

(3) MILES’ (1961) stability condition: A sufficient condition for stability is Ri > 
1/4 everywhere in the flow. 

(4) HOWARD’S (1961) semicircle theorem: The complex wave velocity c of an 
unstable wave with c, > 0 lies in the semicircle 


[c r — (^ 0 max + nomin)/2] 4 C j ^ [(^Omax tiomin)/2] . ($1) 

There are characteristic streamline patterns of fluid motions near critical levels 
known as cat’s eyes. For N = 0 and u 0 — c^z — z c , Kelvin derived from (29) the cat’s 
eye pattern shown in Figure 7a as seen by an observer moving with the wave velocity. 
The existence of closed streamline patterns is a purely kinematic consequence of a fi- 
nite vertical velocity component at the critical level. In a stratified flow (N > 0), the 
cat’s eye pattern becomes asymmetric as indicated by Taylor’s cockeyed cat’s eye in a 
three-layer model fluid (Figure 7b). According to HOWARD and MASLOWE (1973), 
also intermediate forms between Kelvin’s and Taylor’s cat’s eye patterns are possible. 



Fig. 7. (a) Kelvin’s cat’s eye streamline pattern, (b) 

Taylor’s cockeyed cat’s eye streamline pattern. 

There is evidence that cat’s eye structures occur in the time-height distributions of 
MST radar echo power both in the troposphere (ROTTGER and SCHMIDT, 1979) 
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and in the mesosphere (REID et al., 1987) (Figures 8 and 9). The enhanced echo 
power is probably due to enhanced turbulence in thin layers surrounding the regions 
of closed streamlines (KELLY and MASLOWE, 1970). 



Fig. 8. Contour plot of constant MST radar echo power in the time-height plane 
indicating cat's eye structures (from ROTTGER and SCHMIDT, 1979). 
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Fig. 9. Time-height intensity plot of MST radar echo power in three antenna beam 
directions. A cat’s eye is visible between 11:55 and 12:10 LT (from REID et al., 
1987). 


Detailed investigations of cat's eye patterns observed by MST radars have not 
yet been performed, but there are three examples of Kelvin-Helmholtz instabilities 
that were analyzed and compared to numerical model computations (VANZANDT et 
al., 1979; KLOSTERMEYER and RUSTER, 1980, 1981; RUSTER and KLOSTER- 
MEYER, 1983). Figure 10 shows the radial velocity oscillations associated with a 
Kelvin-Helmholtz instability at the bottom side of a polar Jetstream and measured 
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with height and time resolutions of 150 m and 10 s, respectively. The oscillations re- 
veal a period of about 4 min and amplitudes of the order of 1 m s" 1 . The critical level 



Fig. 10. Band pass-filtered time series of radial velocity 
oscillations associated with a Kelvin- Helmholtz instability 
and measured by the SOUSY VHF Radar on 11 April 1980 
(from KLOSTERMEYER and RUSTER, 1981). 



Fig. 11. Measured and computed height profiles of Kelvin- 
Helmholtz-associated velocity oscillations (from KLOSTER- 
MEYER and RUSTER, 1980). 
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is near z = 5.55 km and is indicated by a vanishing amplitude and a sudden phase 
jump of about 100°. Figure 11 shows a comparison between measured and computed 
height profiles of the amplitude and phase of the radial velocity oscillations. Fur- 
ther studies indicate that the growth of the observed Kelvin-Helmholtz instability is 
limited by the onset of static instability producing periodic radar echo power bursts 
(KLOSTERMEYER and RUSTER, 1981). 

4. WAVE- WAVE INTERACTION 

Here we will investigate the general nature of weak resonant interaction between 
wave triplets and will always keep in mind that this is an oversimplification of nat- 
ural processes which consist of resonant and forced energy transfer over continuous 
wavenumber and frequency domains. The actual wave spectrum thus is the result of 
a balance between the effects of sources, dissipation and redistribution due to non- 
linear interaction. Moreover, the application of a weakly nonlinear theory which is 
based on perturbation methods, is restricted to small- amplitude waves. It should be 
mentioned, however, that the stability of a monochromatic internal gravity wave with 
arbitrary amplitude is closely related to weak resonant interactions. This so-called 
parametric instability provides a simple example for the instability of time dependent 
and spatially varying flows and will therefore be discussed in the second part of this 
chapter. 

4.1. WEAK WAVE- WAVE INTERACTION 

If we take into account the nonlinear terms that have been neglected in the previous 
chapters, we obtain a wave equation of the form 

C{w) + Q(w, w) = 0 (52) 

where C and Q are differential operators which are linear and quadratic in uj, re- 
spectively, and w represents any of the perturbation variables. We assume that there 
is a small nondimensional quantity £ characterizing the relative magnitude of the 
nonlinear term in (52). In a straight-forward perturbation procedure we expand w 
in powers of e, 


W — EW\ -h E^W 2 ~b ... (53) 

and solve the equations 

C{w x ) = 0, (54) 

C(w 2 ) = etc. (55) 


Provided the expansion (53) converges we thus obtain a solution for (53) for waves 
of finite but small amplitudes (for details of perturbation methods cf. KEVORKIAN 
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and COLE (1981)). According to Chapter 2, (54) yields plane wave solutions in 
uniform atmospheric models so that 

W] — iW|C lk r_wt 4- complex conjugate (56) 

where the wavenumber vector k and angular frequency satisfy a dispersion relation 


P(k,u/) = 0. (57) 

To second order in e we get from (55) a nonhomogeneous linear equation for w 2 with 
a forcing term depending on W \ . The process can be continued to any order of e. 
Now we consider the case that w x consists of two plane waves, 


W 1 = W a + XVb (58) 

both satisfying (54) and (57). Then (55) becomes 

L(w c ) C(a^a 4" (59) 

where w 2 = u> c is the result of the interaction between w a and w b . Since the quadratic 
operator on the right of (59) produces plane waves with wavenumber vectors k c = 
0, ±2k„, ±2k t , ±(k„ ± k 4 ) and corresponding angular frequencies w c ,u> c has the 
form 


i ^ ^(kc r-wcO q. complex conjugate 
Wc = 2^ Z>(k c>Wc ) 


(60) 


provided D(k c ,u> c ) ^ 0. The summation is performed over all sign combinations 
resulting from the quadratic forcing term in (59). If a = b and both signs are allowed 
for all wavenumbers and frequencies, all cases can be described by 


k a -{- kb -f k c = 0, u a + u >{, + oj c — 0. (61) 

The forced wave e 2 w c plays only a minor role in (53) as compared to e(w a +Wb) unless 
D(k C} u; c ) = 0. Then resonance occurs so that the amplitude of w c grows in space 
and time. Such a resonant wave is also called secular. It dominates the non-resonant 
waves after the interaction has lasted for some time and can therefore be considered 
as the most significant result. 

It depends on the specific form of the dispersion relation whether three interacting 
waves can satisfy the resonance conditions (61) and (57). Clearly, if the waves are 
nondispersive, i.e. propagate at the same speed, all interactions are resonant. For 
internal gravity waves in a Boussinesq fluid with N = const, the resonance conditions 
may be written 


k a ± kb = k c , cos 0 a ± cos 9 b — cos 9 t 


(62) 



287 


using (23). The loci of resonant triads can be calculated from simple geometrical 
considerations and are shown in Figure 12. It should be mentioned, however, that 
the Boussinesq approximation consists of neglecting terms of first and higher orders 
in the ratio of wavelength to density scale height (Chapter 2) and thus is only appli- 
cable if this ratio is much smaller them e (LONG, 1965). 

YEH and LIU (1985) pointed out that resonant wave interaction may play an 
important role in the evolution of the atmospheric wave spectrum. There are in 
particular three classes of resonant triads that can provide efficient energy tremsfer 
between separate parts of the gravity wave spectrum and have first been identified 
by MCCOMAS and BRETHERTON (1977) as elastic scattering, induced diffusion 
and parametric subharmonic instability (Figure 13). By elastic scattering, an upgoing 




Fig. 12. Resonant interaction diagrams for internal gravity waves in a 
Boussinesq fluid for 9 C = 0°,-30° and -60°. Any point on a branch 
defines a resonant triad satisfying (62). The wavenumber vector with 
the least slope is always the vector sum of the other two. Thick (thin) 
branches indicate unstable (stable) triads (after PHILLIPS, 1969). 



Fig. 13. Characteristic classes of resonant triads: 

(a) elastic scattering, (b) induced diffusion, (c) paramet- 
ric subharmonic instability. 
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wave is scattered into a downgoing one by resonantly interacting with a low frequency, 
nearly vertically propagating wave with almost twice the vertical wavenumber of the 
other two waves. This process makes the atmospheric spectrum vertically symmetric 
if it is not so initially. The process of induced diffusion is responsible for the evolution 
of small-scale waves and consists of two almost identical waves which interact reso- 
nantly with a large-scale low-frequency wave propagating almost vertically. Finally, 
in the process of parametric subharmonic instability, two waves with nearly oppo- 
site large wavenumber vectors and almost identical frequencies resonate with a third 
wave having a much smaller wavenumber and twice the frequency. This mechanism 
provides energy transfer from an energetic large-scale wave to small-scale waves at 
subharmonic frequencies. 

4.2. PARAMETRIC INSTABILITY 

There is a close connection between the interaction of weak internal gravity waves 
and the parametric instability of a monochromatic internal gravity wave of arbitrary 
amplitude (MIED, 1976; KLOSTERMEYER, 1982, 1983). For studying parametric 
instability, it is useful to assume a uniformly stratified Boussinesq model and to 
introduce the stream function which describes the two- dimension al motion of an 
incompressible fluid in the x — z plane by 

u = ($ z , - # r ) (63) 

and satisfies (2). We further introduce the buoyancy 

r> P-P* 

B = — g (64) 

where p+ is some suitable constant reference density used in the Boussinesq approx- 
imation to replace p in the inertia terms (e.g. LEBLOND and MYSAK, 1978). 
Buoyancy and Vaisala-Brunt frequency are related by 

iv 2 = B m . (65) 

Equations (3) and (1) then yield 

V 2 ^ + B x = 'J'rV 2 ^ — 'P J V 2 'I' X (66) 

B t - NH Z = <f> x B z - <b,B z (67) 

We note that (66) results from the vorticity equation that is obtained by talcing the 
curl of (3). The plane wave 


'I' = A cosip, B — —N^lw '.A cosy> 


( 68 ) 
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with p ~ lx + mz — urt is a solution of (66) and (67) if the dispersion relation (23) 
is satisfied. The nonlinear terms are identically zero as a consequence of incompress- 
ibility. 

We thus have the opportunity to investigate the stability of a finite- amplitude 
plane internal gravity wave by expressing the stream function and buoyancy fields els 
sums of the basic state (68) and a perturbation: 

'J' = A cos tp + rp ) B = — N 2 lu>~ 1 A cos p + b. (69) 

We define nondimensional variables by 

(x,z) = fc(x,z), 7 = Nt, A: 2 ^" 1 *, B= kN~*B % M = k 2 (2N)~ l A (70) 

with k = (l 2 + m 2 ) 1 / 2 and introduce a rotated coordinate system (£,77) such that 
the r] axis coincides with the direction of propagation of the basic wave (Figure 14). 


i 



Fig. 14. Wavenumber vector of ba- 
sic internal gravity wave, (0,1), and 
Floquet vector (— a sin /?, a cos /?) in 
£ — rj coordinates (from KLOSTER- 
MEYER, 1982). 

Substituting (69) into (66) and (67) and neglecting terms containing products of the 
perturbation quantities 0 and b then yields 

V 2 i p t + sin0^ -f cos 06,, = 2Msin</?(^ + V 2 t/^) ^ 

b t — sin0^ “ cos 6^ = 2 M sin + b$) 

where the tilde of the nondimensional variables has been omitted and V 2 = d 2 /d( 2 + 
d 2 /dr) 2 ,p = T) — cos Ot. The terms on the right of (71) couple the perturbations with 
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the basic state via a periodic coefficient 2M sin so that rj> and b take a parametric 
form dictated by Floquet theory (MIED, 1976): 

(iM) = e , ^ a " n ^ +QCOa/J,7,+Af £ (0„6i)e^. (72) 

js-oo 

As indicated in Figure 14, |a| can be considered as the magnitude of a Floquet vector 
forming an angle j3 with the basic wavenumber vector. Substitution of (72) into (71) 
yields a linear eigenvalue problem for the complex quantities A and (V>;, fy)> 


p\ 'Vj-i + 


(pf =0 

+ r{°V, + (s' 01 - A )bj + rfVi+i + sf ’6j +1 = 0 


(73) 


with j = 0,±1,±2... and constant coefficients p,g,r, s depending on (cf. 

KLOSTERMEYER, 1982). Solutions can be obtained numerically for any given accu- 
racy if the interation equations (73) are truncated at j = ±J with J sufficiently large. 


The numerical search for growing disturbances (Re A > 0) is greatly facilitated by 
considering the case M — ► 0. Then (73) yields for nontrivial solutions (t/> ; , 6 ; ) 


Re A = 0, j cos 0 — ImA ~ 


a cos (9 /?) + j cos 9 

[(a sin 0) 2 + (a cos /? 4- j ) 2 p/2 


(74) 


From (72) and Figure 14, the nondimensional perturbation frequencies and wavenum- 
bers in x — z coordinates are 


— j cos 9 — ImA 

= (a cos(0 -f /?) + j cos 9, asin(0 4- /?) + jsin9) 

so that with (74) 


(75) 


u)j — cos 9j (76) 

where 9 } is the angle between and the horizontal. For M — + 0, the disturbances 
thus are freely propagating waves that satisfy the resonance conditions 

k;+i - k, = k, u j+ 1 - Uj = w (77) 

for aJl j. Computations for M > 0 show in agreement with HASSELM ANN’s (1967) 
criterion for nonlinear wave stability that ReA > 0 (i.e. instability) for resonant sum 
interactions and ReA = 0 (i.e. stability) for resonant difference interactions. In Fig- 
ure 12, the loci of unstable triads are indicated by the thick branches. 

In a laboratory experiment, DAVIS and ACRIVOS (1967) demonstrated that an 
internal gravity wave propagating along a diffuse stratified interface between fluids of 
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different densities may become distorted by the growth of a three- wave resonance giv- 
ing rise to local turbulent mixing. Their photographs are shown in Figure 15. MCE- 
WAN and ROBINSON (1975) found good agreement between predicted and observed 
small-scale parametric instabilities in a cylindrical chamber. Although the potential 







Fig. 15. Progressive distortion of an internal gravity 
wave by resonant wave- wave interaction (from DAVIS 
and ACRIVOS, 1967). 


role of parametric instability in forming internal wave spectra in the atmosphere has 
been stressed by some authors (MCEWAN and ROBINSON, 1975; DRAZIN, 1977), 
there are only few atmospheric observations that have been interpreted in terms of res- 
onant triads or parametric instability. Thus the power spectrum of a strong internal 
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Fig. 16. Doppler shift of a 8 MHz continuous radio wave 
recorded at four Japanese stations after the eruption of 
Mount St. Helens on 18 May 1980.' In the lower three 
records, an internal gravity wave with a period near the 
local Vasisala-Brunt period is strongly disturbed by short- 
period oscillations after 22:55 UT (from LIU et ah, 1982). 

gravity wave in the lower thermosphere which was generated by am eruption of Mount 
St. Helens, was explained in terms of parametric instability (KLOSTERMEYER, 
1984). Figure 16 showing the Doppler shift of an 8 MHz radio wave reflected from 
the ionosphere near 200 km height, indicates strong high-frequency disturbances su- 
perposed on a gravity wave at Kyoto, Yokaichi and Uji after 2255 UT. Note a certain 
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resemblence of these records to Figure 15c. 

ROTTGER (1987) has noted that time series of spectra intensity plots obtained 
from MST radar measurements also show high-frequency oscillations superposed on 
low-frequency gravity waves. An example is seen in Figure 17 at 69.6 km after 
1300 AST. Rottger further points out that the amplitude the low-frequency waves 
does not grow with inreasing height (visible between 1200 and 1300 AST in the height 



1200 1300 1400 1600 AST 


Fig. 17. Time series of spectra intensity plots measured with the mo- 
bile SOUSY VHF Radar at the Arecibo Observatory (from ROTTGER, 

1987). 

interval from 69.6 to 76.8 km). Since wave breaking discernible from strong echo in- 
tensities does not occur some other saturation process such as parametric instability 
must be considered. 

Short- period gravity waves with frequencies close to the local Vaisala-Brunt fre- 
quency occur frequently in mesospheric MST radar observations of vertical or near 
vertical radial velocities (HARPER and WOODMAN, 1976; MILLER et al., 1978). 
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These waves do not show any vertical phase variation over height intervals of sev- 
eral kilometers indicating that they propagate horizontally and are generated in situ 
rather than in lower or higher atmospheric regions. Occasionally the local source 
mechanism is Kelvin-Helmholtz instability of a wind shear due to long-period iner- 
tial or tidal waves (KLOSTERMEYER and RUSTER, 1984; YAMAMOTO et al., 
1988). But in general the simultaneously observed wind shear is too weak for Kelvin- 
Helmholtz instability to set in. A numerical simulation of the wave motions shown in 
Figure 3 indicates that parametric instability could be another source candidate of 
the observed short-period waves. Figure 18a shows the same observations as Figure 
3 but with very high-frequency noise and high-frequency waves removed by recursive 
Butterworth filters (continuous and dashed curves respectively). Figure 18b shows 
the numerical simulation where the dashed lines represent the finite-amplitude basic 
wave and the continuous lines the sum of the basic wave and the fastest growing para- 
metric instability mode. Note that the height interval in Figure 18b is 75% of the basic 
vertical wavelength whereas in Figure 18a, it is only about 10%. The observed and 
computed short-period oscillations show good qualitative agreement. Both reveal in 
particular considerable temporal amplitude modulation, and the dominating period 
is not an integer multiple of the basic period. There is in general also no vertical phase 
variation. But the computed short-period oscillations show sudden phase reversals, 
e.g. at Nt = 20 and kz = 4.8. The loci of sudden phase reversals lie on basic wave 


SOUSY VHF RADAR 
21 JAN 19BA 



Fig. 18a. Same as Figure 3 but with noise and 
high-frequency waves removed (continuous and dashed 
curves, respectively). The time series at z = 70.6 km 
is not continued beyond 12 LT with regard to very 
small signal-to-noise ratios. 
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Nt 

Fig. 18b. Numerical simulation of a paramet- 
ric instability. The dashed and continuous lines 
represent the basic wave and the sum of basic 
wave and fastest growing instability mode, re- 
spectively. Time, height and vertical velocity 
component are normalized according to (70). 


fronts. Whereas there is no clear evidence for phase reversals in Figure 18a (perhaps 
due to the small height interval ) they become visible if velocity measurements can 
be obtained over larger height intervals (e.g. Figure 1 of MILLER et al. (1978)). 
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